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Abstract 

The global structure of solutions of the Einstein equations cou¬ 
pled to the Vlasov equation is investigated in the presence of a two- 
dimensional symmetry group. It is shown that there exist global CMC 
and areal time foliations. The proof is based on long-time existence the¬ 
orems for the partial differential equations resulting from the Einstein- 
Vlasov system when conformal or areal coordinates are introduced. 


1 Introduction 

In general relativity the gravitational field is expressed in terms of the geom¬ 
etry of spacetime. The matter constituting self-gravitating physical systems 
is described using certain matter fields which are also geometrical objects. 
The spacetime metric is required to satisfy the Einstein equations and these 
are coupled to equations of motion for the matter fields. The resulting 
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Einstein-matter system of partial differential equations is the central math¬ 
ematical element of the theory. Thus to understand the mathematical con¬ 
tent of general relativity we require an overview of the solutions of these 
equations and their qualitative behaviour. 

When different kinds of physical situations are considered solutions of 
the Einstein-matter equations with certain boundary conditions or spatial 
asymptotics will be particularly relevant. One choice which avoids the issue 
of boundary conditions or asymptotic conditions in the strict sense is to give 
initial data on a compact manifold. This possibility, which is that studied 
in the present paper, is of interest for applications to cosmology. In practice 
it leads to imposing periodic boundary conditions. 

A further choice which has to be made in order to get a concrete math¬ 
ematical problem is to decide which kind of matter fields to consider. Here 
we choose collisionless matter, where the equation of motion is the Vlasov 
equation. The special interest of this kind of matter has been discussed 
in several places (see |R13|, | R15(| , []A2(] ) and the relevant facts will not be 
repeated here. A fundamental insight is that collisionless matter acts as a 
source of the gravitational field in a way which allows specihc features of the 
matter model to remain in the background while revealing basic properties 
of the dynamics of self-gravitating matter. A fact on the PDE level related 
to this is that in many cases long-time existence theorems for solutions of 
the Einstein-Vlasov system are to be expected. For other matter models 
coupled to the Einstein equations this is not true. Results on formation of 
singularities in finite time in the case of one matter model, dust, can be 
found in H and [p. 

While the Einstein-matter system itself is naturally expressed in geomet¬ 
rical terms the application of the theory of partial differential equations to 
study it requires the geometry to be parametrized in a suitable way by the 
use of coordinates or other auxiliary constructs. An important step towards 
investigating the dynamics of solutions is to find parametrizations which 
are at once practical and applicable in sufficient generality. Of particular 
significance is finding a function which can act as a good time coordinate in 
the situation of interest. This task is the focus of the following. It will be 
shown that under appropriate symmetry assumptions it can be solved in a 
very satisfactory way. 

A central role in solving the geometric problems addressed in this paper 
is played by two long-time existence results for certain systems of PDE in one 
space dimension. One of these is a global in time existence theorem while 
the other is a continuation criterion which says that solutions continue to 
exist as long as a certain quantity does not vanish. A number of techniques 
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used in the proofs are adapted from known arguments for the Vlasov-Poisson 
and Vlasov-Maxwell systems while others are specific to the Einstein-Vlasov 


case. 


Consider a solution of the Einstein-Vlasov system evolving from initial 
data on a compact spacelike hypersurface. This paper is mainly concerned 
with the case where the initial hypersurface is a three-dimensional torus with 
periodic coordinates {6,x,y). Moreover solutions are considered where both 
the metric and the matter fields are invariant under translations in the x and 
y directions. Taking account of the periodic identifications involved we see 
that these solutions admit a two-dimensional symmetry group isomorphic to 
a two-torus x S^. These spacetimes will be referred to as sym¬ 

metric. No additional symmetry assumptions will be imposed. Analogues of 
the results of this paper for various cases with higher symmetry, including 
certain subcases of symmetry, have been obtained previously. A survey of 
these earlier results can be found in the introduction of HARR |. The present 
paper can be seen as the culmination of a development concerning global 
geometrically defined time coordinates in solutions of the Einstein-Vlasov 
system with at least two symmetries. 

On a symmetric spacetime it is possible to define a function R, the 
area function, as follows. If p is a point of the spacetime then R{p) is equal 
to the area of the orbit of the action of which contains p. Evidently 
the function R is itself invariant under the action of . We say that a 
symmetric initial data set is flat if the spacetime gradient of i? in a Cauchy 
development of this data set vanishes identically on the initial hypersurface. 
Otherwise it will be called non-flat. Whether an initial data set is flat 
in this sense can be determined intrinsically from the initial data without 
having to know anything about a Cauchy development. The reason for the 
terminology is that, as follows from [ [Rll[ |, if an initial data set is flat in this 
sense any Cauchy development of it is flat in the sense that its Riemann 
curvature tensor vanishes everywhere. Cf. the discussion in Section 

The main theorems will now be stated. The first concerns the existence of 
a global time coordinate of constant mean curvature. Consider a spacetime 
which evolves from initial data on a compact hypersurface. In the following 
it will always be assumed implicitly that the initial datum for the particle 
density has compact support. A real-valued function t on the spacetime is 
called a constant mean curvature (CMC) time coordinate if each of its level 
hypersurfaces is compact and spacelike, it has constant mean curvature and 
the value of the mean curvature there is equal to t. 


Theorem 1 Let f) be the maximal globally hyperbolic development 
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of non-flat C°° initial data for the Einstein-Vlasov system with symme¬ 
try. Then M can he covered by compact spacelike hypersurfaces of constant 
mean curvature with each value in the range (— 00 , 0 ) occurring as the mean 
curvature of precisely one of these hypersurfaces. 

The second theorem concerns the existence of an areal time coordinate. 
Consider a spacetime with symmetry. A real-valued function t on the 
spacetime is called an areal time coordinate if each of its level hypersurfaces 
is compact and spacelike and the value of t on the hypersurface is everywhere 
equal to that of the area function R. 

Theorem 2 Let {M, g^i^, f) be the maximal globally hyperbolic development 
of non-flat C°° initial data for the Einstein- Vlasov system with symme¬ 
try. Then M can he covered by compact spacelike hypersurfaces of constant 
area function R with each value in the range {Ro,oo) occurring as the value 
of the area function on precisely one of these hyper surf aces. Here Rq is a 
non-negative real number. 


Notice that these theorems include the vacuum case / = 0. Theorem 2 
was proved in the vacuum case in [ BCIM |. Under the additional assump¬ 
tion of Gowdy symmetry, which consists in augmenting the action by 
a suitable reflection symmetry. Theorem 2 was proved in |^|. In | ARR ] 
an argument was sketched which indicates that Theorem 1 also holds under 
the assumption of Gowdy symmetry. The results of both theorems extend 
to the case of local U{1) x U{1) symmetry, as defined in Rll |. Since there 
is no essential difference in the proofs this generalization will not be men¬ 
tioned further. The fact that the theorems are stated for (7°° initial data is 
a matter of convenience. Keeping track of derivatives would allow analogous 
results to be proved for initial data of finite differentiability. 

The paper is structured as follows. The next section contains some 
definitions which will be needed. The two sections after that contain the 
basic PDE analysis in the contracting and expanding directions respectively. 
The proofs of the main theorems are given in Section 


2 The Einstein-Vlasov system with symmetry 

Gonsider the manifold M = M x T^. Let act on in the obvious way, 
arising from the action of on with Gartesian coordinates {6,x,y) by 
translations in x and y. Gorrespondingly T^ acts on M with coordinates 
{t,9,x,y). A spacetime with underlying manifold M defined by a metric 
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ga 0 and matter fields is said to be symmetric if the metric and matter 
fields are invariant under the action of T^. The orbits of the group action 
will be referred to as surfaces of symmetry and a hypersurface will be called 
symmetric if it is a union of surfaces of symmetry. It is now clear how 
to define abstract Cauchy data for the Einstein-matter equations with 
symmetry. They should be defined on and invariant under the action of 
T^. It will be assumed throughout that the metric and the matter fields are 
C°°. 


Suppose now that a matter model is chosen for which the Cauchy prob¬ 
lem for the Einstein equations is well-posed. The example of interest in 
the following is that of collisionless matter satisfying the Vlasov equation. 
Corresponding to initial data for the Einstein-matter equations with sym¬ 
metry there is a maximal Cauchy development. We will construct a certain 
local coordinate system on a neighbourhood of the initial hypersurface in the 
maximal Cauchy development. On the initial hypersurface itself we choose 
periodic coordinates { 6 ,x,y) as above. The isometries are given by transla¬ 
tions in X and y. These coordinates can be extended uniquely to a Gaussian 
coordinate system {t, 9, x, y) on a neighbourhood of the initial hypersurface. 
On general grounds the action of on the initial data extends uniquely to 
an action on the maximal Cauchy development M by symmetries (see e.g. 
[FR|, Section 5.6). Gauss coordinates inherit the symmetries of their initial 
values and the spacetime. Hence the components of the metric in these co¬ 
ordinates depend only on t and 9. It is then a matter of simple algebra to 
see that in these coordinates the metric can be written in the form 


- dt'^ + d9^ + e^^ [dx + Ady + {G+AH)d9f + [dy + Hd9] 2.(1) 

for functions (fj, U, A, G, H, R) of t and 9 which are periodic in 9. It can also 
be read off that the function R in this form of the metric coincides with the 
area function mentioned in the introduction. 

A region covered by coordinates of this type which is invariant under the 
action of can be quotiented by the group action to get a two dimensional 
quotient manifold Q coordinatized by t and 9. The quotient inherits a 
Lorentzian metric by the rule that the inner product of two vectors on the 
quotient is equal to the inner product of the unique vectors on spacetime 
which project onto them orthogonally to the orbit. On Q we can pass to 
double null coordinates (u, v) on a neighbourhood of the quotient of the 
initial hypersurface, i.e. to coordinates whose level curves are null. Defining 
new coordinates hy t = ^{u — v) and 9 = ^{u + v) puts the metric on Q 
into conformally flat form. By pull-back these define new coordinates on M 
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where the metric takes the form 

(2) 

for functions {rj, U, A, G, H, R) of t and 9. A coordinate system of this type is 
called a conformal coordinate system. It has now been seen that conformal 
coordinates always exist on some neighbourhood of the initial hypersurface 
in a spacetime evolving from data with T‘^ symmetry. In addition it is 
possible to choose the double null coordinates in such a way that the initial 
hypersurface coincides with t = 0 and the metric coefficients are periodic in 
9. 

To conclude this section we formulate the Einstein-Vlasov system which 
governs the time evolution of a self-gravitating collisionless gas in the context 
of general relativity; for the moment we do not assume any symmetry of the 
spacetime. All the particles in the gas are assumed to have the same rest 
mass, normalized to unity, and to move forward in time so that their number 
density / is a non-negative function supported on the mass shell 

PM := {'q^yV^v'' = -1, > O} , 

a submanifold of the tangent bundle TM of the space-time manifold M with 
metric gap- Here denotes the components of the Minkowski metric and 
denote the components of a tangent vector in an orthonormal frame e^. 
We use coordinates (t, x“) with zero shift and the frame components u* to 
parametrize the mass shell; Greek indices always run from 0 to 3, and Latin 
ones from 1 to 3. On the mass shell PM the variable becomes a function 
of the u*: 

The Einstein-Vlasov system now reads 

^'%(/) - = 0 

„ ,, ^ dv^dv'^dv^ 

v^v f - 

-Vo 

where are the Ricci rotation coefficients, the Einstein tensor, and 
the energy-momentum tensor. In the case of symmetry the number 
density / is assumed to be invariant under the action induced on PM by 
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the action of on M. The orthonormal frame used to parametrize PM in 
the following is also assumed to be invariant. As a result / is independent of 
X and y and is a function of the variables (t, 9, u^). An explicit choice 

of invariant orthonormal frame for the metric (^) is given by 



3 Analysis in the contracting direction 

As we will see in more detail later, a non-flat symmetric solution of 
the Einstein-Vlasov system represents a cosmological model which, after a 
suitable choice of time orientation, evolves from an initial singularity and has 
a phase of unlimited expansion at late times. In this section we consider the 
Einstein-Vlasov system in the contracting direction (i.e. evolving towards 
the initial singularity) and use conformal coordinates in which the metric 
takes the form (^). That coordinates of this type can always be found 
near the initial hypersurface was shown in the last section. The invariant 
orthonormal frame adapted to these coordinates exhibited there will be used 
in this section. The explicit form of the Einstein-Vlasov system will now be 
given. We set T = Gt + AHt- The quantities T and Ht will be referred to 
as twist quantities. 

The Einstein-matter constraint equations 



( 3 ) 


( 4 ) 


The Einstein-matter evolution equations 



Q-2ri+W 


2 


( 5 ) 
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Rtt — Ree 
m - Vee 


- Pi) + 

p4f/ 

= 

3i?2e-2’? 


j^^- 2 r,+iU 


-r^ + 


i?3e-2»? 


-m 


t ’ 


(4 - Al) - 


^- 2 r)+iU 


-H? - 


(6) 

( 7 ) 

( 8 ) 


Auxiliary equations 


de[Re-^^+^^T] = -2Re^J2, (9) 

dt[Re-^^+^^T] = 2Re'^Si2, ( 10 ) 

deiR^e-^^Ht) + AgT = -2R^e^-^^ J 3 , (11) 

dt{R^e-^^Ht) + Re-^^+^^AtT = 2 R^e'^-^^ S 13 . ( 12 ) 


The Vlasov equation 


df df 
'di ^ 


ive - Ue)v^ + {pt - Ut)v^ - 


+(t^» - f )^ - + RH,v^) 


., 2„3 


dv^ 


1„.2 


Utv^ + Ue- 


V V 


K 

(9^2 


M. 


Rg . 


V,3 


(^ - Ut)v^ - {Ug - ^) 


R 

The matter quantities 

p{t,e) = 
Pk{t,e) = 
Jk{t,0) = 


R' 


+ 




R 


■{At + Ag—^) 


K 

dv^ 


v^f{t, 6 , v) dv 


n 


k \2 


/ k 3 V 

r 

,,fc 


-f{t,9,v)dv, A: = 1,2, 3 


/ 

Jm 


v^f{t, 6 , v) dv 

,0 f{t,0 ,v)dv 


= 0. (13) 

(14) 

(15) 

(16) 


yiy^ 


K3 V 


Sjk{t, 6 ) 


(17) 



Let a smooth symmetric solution of the Einstein-Vlasov system writ¬ 
ten in conformal coordinates be given on some time interval We 

want to show that if this interval is bounded and if R is bounded away from 
zero there then /, R, r], U, A, G, H and all their derivatives are bounded 
as well, with bounds depending on the data at t = to the lower bound 
on R, and that the supremum of the support of momenta at time t, 

Q{t) := sup{|u| : 3(s,0) e [t,to] x 5'^such thatf{s,9,v) / 0}, (18) 


is uniformly bounded. Note that these conditions imply that the matter 
quantities and their derivatives are uniformly bounded. In proving these 
statements it is assumed that the initial data at t = to are non-flat, as in 
the assumptions of the main theorems. 

The description of the proofs in this section is modelled on that of [Al] 
and highlights the places where there are differences due to non-vanishing 
twist. 

Step 1. (Monotonicity of R and bounds on its first derivatives.) 

This is a key step and follows the arguments in [BCIM]. We have to check 
that the matter terms have the right signs so that these arguments still hold. 
The bounds on R and its first derivatives will play a crucial role when we 
control the matter terms below. 

First we note that Vi? is timelike. This is a consequence of Proposition 
3.1 of [nil]. A symmetric solution of the Einstein-Vlasov system satisfies 
the hypotheses of that proposition. Since the initial data set is non-flat the 
Hawking mass is non-zero somewhere. Then the proposition implies that it 
is non-zero everywhere. As a consequence Vi? is everywhere timelike. 

Next we show that dtR and \d 0 R\ are bounded into the past. Let us 
introduce the null vector fields 


1 


di = -^{dt + de) 




(19) 


and let us set Eg = F\ = dxF for any function F. The evolution 

equation (^) can be written 

dxR^ = - Pi) + + R^Hf), (20) 

or equivalently, 

d^Rx = |e2(’'-^)(p - Pi) + :?^(e4C^r2 + R^H^). (21) 
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The right hand side is positive since p > Pi and we can conclude, arguing 


as in Step 1 in Section 4 of |M|, that both Rt and [ii^l are bounded into the 
past. Hence R is uniformly bounded to the past of the initial surface. 

Step 2. (Bounds on U, A and p and their first derivatives.) 

The bounds on Ut-,At,Ug and Ag to the past of the initial surface are ob¬ 
tained by a light-cone estimate, which in this case, with one spatial dimen¬ 
sion, is an application of the Gronwall method on two independent null 
paths. Then, by combining these results, one obtains the desired estimate. 

The functions involved in the light-cone argument are quadratic func¬ 
tions in the first order derivatives of U and H, defined by 


y = RUtUe + —AtAe. 


( 22 ) 

(23) 


We will see below that if we let the vector fields along the null paths act on 
X + Y and X — Y we obtain equations appropriate for applying a Gronwall 
argument. 

Let us now derive bounds on U and A and their first order derivatives. 
By using the evolution equations (I) and (|6|) we find 


dx{X + Y) = 


-1 

2 V 2 

R 


,4(7 


Re Ut - Ui + 


AR? 


{-Aj + AD 


+ e2(^-^)(p - Pi + P2 - P3)) 

pPU 

+—H5(p2e2(^-'7)rPi + 2Pe2(’?-f^)523), 

ZK 


and 


d^{X-Y) = 


-1 

2V2 

R 


p4(7 

R^[U^-U^ + —i-Al + Al) 


+^Ux + e2(^-^)(p - Pi + P2 - P3)) 


p2(7 

+—Ax{R^e^^^-^^THt + 2 Re^^^-^'^ S23). 

IH 

It turns out that X and Y can be bounded by integrating these equations 
along null paths starting at a general point {ti,6) in the past of the initial 
hypersurface and ending at the initial to surface t = to. The boundedness 
of the integrals which arise follows from the equations (^) and (21) and 
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the boundedness of and Rx in a similar way to the corresponding step 
in Section 4 of ||A1|| . Up to multiplication by functions which are already 
bounded and functions which can be bounded linearly in terms of X the 
terms involving the twist whose integrals must be estimated are 
and The first of these, multiplied by i?/4, occurs as one of 

the positive summands on the right hand side of (|^ ) and (^). Since by 
assumption R is bounded below the first term can be controlled. For the 
second we can use the elementary inequality 




(24) 


and the occurrence of 9'S a summand in (^) and (^). Bounds 

on X and Y follow by applying Gronwall’s inequality as in [Al]. Thus, as 
long as R stays uniformly bounded away from zero we conclude that U and 
its first order derivatives, and thus also A and its first order derivatives, are 
bounded. Bounds on |7/|, |? 7 t| and \riQ\ are obtained in a similar way since the 
evolution equation (P) can be written 


,4(7 


28 x 11 ^ = Ui - Ut + -^{Ai - Ai) - 


2?j+4(7 


or equivalently. 


,4(7 




3/?2e-2’? 


,-2r?+4(7 


2d^r^x = U^B-Ut + ^{At-Ai)- 


-r^ - 


3i?2e-2^ 


e2G-^)p3, 

(25) 

H? - 


(26) 

The integrals of the right hand sides of these equations along null paths are 
bounded since P 3 < p — Pi- The terms involving T and Ht can be handled 
as above. Thus we hnd that rj is uniformly bounded to the past of the 
initial surface as long as R stays bounded away from zero. 

Step 3. (Bound on the support of the momentum.) 

Note that a solution / to the Vlasov equation is given by 


/(i, v) = /o(0(o, t, e, v), V (0, t, 6, v)), 
where 0 and V are solutions to the characteristic system 

('T/ 2'\2 

= -{m-Ue)V^-{pt-Ut)V^ + Ue^-^ 


(27) 


11 
















dV'^ 

ds 

dV^ 

ds 


ReAV^r , Ae 

-(f'o --R )—+-R <= — 


= -UtV^-Ue- 


FO 


^R '* ® R^ VO 


^2U 




and Q{s,t, x,v), V{s,t, x,v) is the solution that goes through the point 
{6, v) at time t. Let us recall the definition of 

Q{t) := sup{|u| : 3(s,0) G [t, to] x <S'^such that/(s,0,u) / 0}. 

We also define , j = 1,2,3 in the obvious way where |u| is replaced by 
|u-^ |. If Q{t) can be controlled we obtain immediately from (^)- dT^) bounds 
on p,Jk,Pk, and Sjk, j,k = 1,2,3, j / k, since ||/||oo < ||/o||oo from (p7|). 
First we note that the quantities 


e^F^ Ae^V'^ + Re-^V^ 

are conserved which a simple computation shows by using the equations for 
dV'^/ds and dV^/ds above. More generally this is a consequence of the fact 
that if 7 is a geodesic and k a Killing field then k) is conserved along 
the geodesic. Here 'y' is the tangent vector to 7 and since the particles follow 
the geodesics of spacetime the tangent vector can be expressed in terms of 
v^. We chose k = dx and k = dy as Killing fields to derive the conserved 
quantities given above. Since U and A are uniformly bounded as long as 
R stays bounded away from zero we conclude that and and thus 
and are bounded as well. In order to bound Q{t) we need to control the 
remaining component Q^. From the auxiliary equations (|l^) and (jl^) we 
conclude that F and Ht can be bounded by the quantities involving tj, U, At, 
which are known to be bounded, and by | 5 i 2 | and jS'isj. Since and are 
controlled we immediately get from ( [T7| ) that | 5 i 2 | and | 5 i 3 | are bounded 
by a constant times Q^. We conclude that as long as R stays bounded away 
from zero the terms involving F and Ht in the characteristic equation for 
dV^/ds can be estimated by C{t)Q^{t), where C{t) is uniformly bounded 
on closed time intervals. Now, since the field components U, A and 7 and 
their first derivatives are known to be bounded on (t-, to] (as long as R stays 
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bounded away from zero) we obtain from the characteristic equation for 

Q\s)ds<\Q\to)\+C{t) Q\s)ds, t<to. 

Jt Jt 

Note that (5^(to) is bounded by a positive constant since /o has compact 
support. This inequality leads to a Gronwall inequality for Q^{t) and we 
conclude that Q^{t) is uniformly bounded on [t, to]- 

Thus all the field components, their first derivatives and the matter terms 
are known to be bounded on (t_, to]) as long as R stays bounded away from 
zero. 

Step 4- (Bounds on the second order derivatives of the field components 
and on the first order derivatives of /.) From the Einstein-matter constraint 
equations in conformal coordinates we can express RtQ and Rqo in terms 
of uniformly bounded quantities, as long as R stays bounded away from 
zero. Therefore these functions are uniformly bounded and equation ( 0 ) 
then implies that Ru is uniformly bounded as well. 

In the vacuum case one can take the derivative of the evolution equa¬ 
tions and repeat the argument in Step 2 to obtain bounds on second order 
derivatives of U and A. Here we need another argument. First we write the 
evolution equations for U and A in the forms 


Utt-Uee = + 


2R 

^4U 


2R 


p-2r}+4U 1 

+ ^(^t - Ae){A, + Ae) + r" + 


and 


Att — Agg = 


ee 


iAe + A,) + ^ 

-2(A, - Ai,)(U„ + U,) - 2(A, + A,)(U, - Ue) 

+R^e-^^THt + 2i?e2('?-2^)S23, 


where k denotes p — Pi+ P 2 — P 3 . Taking the 0-derivative of these equations 
gives 


dxd(Ug 


T) D p4[/ 

+ A^dxAg) 

+ (28) 
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and 


dxd^Ag — L + ~ “^U^dxAg + 2Axd^U0 + 2Uxd^AQ 

+2A^dxUe + R^e-^\THt)e + 2Re^^^-^^') { 823 ) 9 , (29) 

Here, L contains only T,Ht,K and S 23 , first order derivatives of U,A and 
r], and first and second order derivatives of R, which all are known to be 
bounded. These equations can of course also be written in a form where the 
left hand sides read d^dxUg and d^dxAg, respectively. By integrating these 
equations along null paths to the past of the initial surface, we get from a 
Gronwall argument a bound on 


sup(|(9^C/0| + \dxUe\ + \d^A9\ + |(9a-40|), 
eesi 

as long as R is bounded away from zero, under the hypothesis that the in¬ 
tegral of the differentiated terms Tg^Htg^ng and { 823)9 can be controlled. 
That the first two terms are bounded follows immediately from the auxiliary 
equations (^) and (plj ) in view of the bound on Q{t). In order to bound the 
matter terms we make use of a device introduced by Glassey and Strauss 
[]GS[] for treating the Vlasov-Maxwell equation. Since there is no essential 
difference to the case treated in |A1| we do not repeat the proof here. This 
procedure allows the integrals of the differentiated matter terms to be con¬ 
trolled and the Gronwall argument referred to above goes through. So we 
obtain uniform bounds on \d^U 9 \,\dxU 9 \,\d^A 9 \, and \dxA 9 \, and therefore 
also on \U 99 \,\Ute\-, and \At 9 \, as long as R is bounded away from zero. 
The evolution equations (P) and (P) then give uniform bounds on \Uu\ and 
I Att I . Bounds on second order derivatives of G and H follow from the auxil¬ 
iary equations. By differentiating equation (^), it is now straightforward to 
obtain bounds on the second order derivatives of rj, using similar arguments 
to those already discussed here, in particular the integrals involving matter 
quantities can be treated as above. Bounds on the first order derivatives of 
the distribution function / may now be obtained from the known bounds 
on the field components from the formula 


f{t, 6 , v) = /o(0(O, t, 0,v),V (0, t, 9, v)), 


(30) 


since /o is smooth and since dQ and dV (here d denotes 81,89 or 8 ^) can 
be controlled by a Gronwall argument in view of the characteristic system 
and the auxiliary equations (|^-(]^). 
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Step 5. (Bounds on higher order derivatives and completion of the proof.) 
It is clear that the method described above can be continued for obtaining 
bounds on higher derivatives as well. Hence, we have uniform bounds on the 
functions R, U, A, rj and / and all their derivatives on the interval (t_, to] if 
R > e > 0. This implies that the solution extends to t —> —oo as long as R 
stays bounded away from zero. 

Later we will require a slight generalization of these results in order to 
show that the arguments of Section 5 of |ARR| generalize to cover the case 
of symmetry. Once it has been established that Vi? is timelike, the 
estimates in the later steps hold for any subset Z of the half-plane t < Iq 
provided Z is a future set. By definition this means that any future directed 
causal curve in the region t < to starting at a point of Z remains in Z. (For 
information on concepts such as this concerning causal structures see e.g. 
[|HEj ].) Thus if R is bounded away from zero on Z and t is bounded on Z 
then all the unknowns and their derivatives can be controlled on Z. 

Now consider a special choice of the subset Z, namely that which is 
defined by the inequalities ti < t < to and 6 i + to — t < 9 < 62 — to +1 for 
some numbers 6 * 1 , O 2 and ti satisfying the inequalities 6*1 < 62 and ti>to~ 
(1/2) (02— 01 ) • Suppose a solution of the Einstein-Vlasov system in conformal 
coordinates defined on this region is such that R is bounded away from zero. 
Then the functions defining the solution extend smoothly to the boundary 
of Z at t = ti. They define smooth Cauchy data for the Einstein-Vlasov 
system. Applying the standard local existence theorem (without symmetry) 
allows the solution to be extended through that boundary. Repeating the 
construction of conformal coordinates in Section ^ then shows that we get 
an extension of the solution written in conformal coordinates through that 
boundary. 


4 Analysis in the expanding direction 

In this section we want to investigate the Einstein-Vlasov system with 
symmetry in the expanding direction. We write the system in areal coordi¬ 
nates, i.e., the coordinates are chosen such that R = t. The circumstances 
under which coordinates of this type exist are discussed in Section p. We 
prove that for initial data on a hypersurface of constant time corresponding 
solutions exist for all future time with respect to the areal time coordinate. 
In order to extend a solution defined on a finite time interval in these coor¬ 
dinates to one which exists globally in time it is sufficient to obtain uniform 
bounds on the field components and the distribution function and all their 
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derivatives on a finite time interval [to, t+) on which the local solution exists. 
For in this case the functions defining the solution are uniformly continuous 
and extend continuously to t = t+. Since all derivatives are bounded the 
extension is smooth and we obtain new initial data on t . Applying the 
general local existence theorem for the Einstein-Vlasov system |CB| allows 
the spacetime to be extended beyond f Conformal coordinates can be 

introduced on a neighbourhood of f and from those areal coordinates 
can be constructed. In this way an extension to the future of the solution 
written in areal coordinates is obtained. 

Below the form of the metric and the Einstein-Vlasov system are given 
in areal coordinates. The functions a, rj, U, A, G, H all depend on t and 6 
and the function / depends on t, 6 and u G and as before we have set 
r = Gi -|- AHt. The orthonormal frame used to parametrize PM is 


a 




dt 


89 dx dy 


,e 


— eM- 
ax 


— -A— 
dy dx 


Metric 


9 


e2(’?-^)(-adt2 + ^02) ^ ^ + AH)def 

+e-^^t\dy + Hdef (31) 


The Einstein-matter constraint equations 


m 

t 

t 

at 


Ut -\- allg -\- 


pW 

-^{Aj + aA^) + 


- 2 »? 


+eHt) 


4U 

2UtUg + ^AtAg - 

2ta‘^e^^^-^^Pi - p)- ate-2^(e^^r2 + 


(32) 

(33) 

(34) 


The Einstein-matter evolution equations 


Ptt — onrjgg = 


pgag Ptat ag 


+ 


agg 


2a 




,4C/ 




- - ie-2^e"^r 
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( 35 ) 




Utt — aUeo = 


Att — aAeg = 


Ut , Ugae , Utat , ^ ,2 „ 42 


eJ 


3-277 4(7 


+ 


-A + 


32(t7-(7) 


a 


(p — A + A — A 


2 2 

, ae^e , 4 4 rr 1 4 4 rr 

— H-7^-h - ^AtUt + AaAgUg 

t I la 

+t^e-‘^^THt + As- 


(36) 


(37) 


Auxiliary equations 


A[e 


-2r,^-l/2^4[/-pi ^ 


r] = 


9i[e-2'?to-V2e4^r] 

A[e"^A-^/ 2 (^g 4 C/p ^ 
9t[e-2^to-^/2(^g4C/p ^ 


-2eV2, (38) 

2ta^/A^5i2, (39) 

-2e^AJ2 - 2te^-^^J3, (40) 

2ta^/‘^e^{ASu + te-^^Sia)- (41) 


The Vlasov equation 

^ y/av^ df 
dt 89 

Aq v'^v^ y/aUg 


+ 


{V0 -Ug + ^)\/an° + (r/i - Ut)v^ 
la 

- (A)2) + e-^(e2^rA + tHtv^) 


dv^ 


Utv"^ + y/aUg 


v^v'^ 


K 

dv‘^ 




( _ Uty - y/^Ug^ + 


{At + y/aAg — ) 


K 

dv^ 


= 0. (42) 


The matter quantities are defined as in (O)-®- 
Step 1. (Bounds on a,U, A,G, H and iy.) 

In this step we first show an “energy” monotonicity lemma and then we 
show how this result leads to bounds on ?y := ?y + In a/2 and on U and A. 
Let E{t) be defined by 


E{t) 


L 


1 e^^ 1 

a ‘^Ut + yfaUg +-^{a + yfaA^) 


+ - 


51 


(e^^r 2 + t^H^) + y/^e?^-^-^)p 


de. 


(43) 
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Lemma 1 E{t) is a monotonically decreasing function in t, and satisfies 


A. 

dt 


m = - 


L 


uf 


si L 
a 


p^U p-2ri 

+ + . j— 

a 4ya 


(e' 


4U 


+ 2t‘^Hf) 


+^e^^^-^\p + P,) 


de < 0 . 


(44) 


Proof. This is a straightforward but lengthy computation. Let us sketch the 
steps involved. After taking the time derivative of the integrand we use the 
evolution equations for U and A to substitute for the seond order derivatives, 
we use the auxiliary equations to express second order derivatives of G and 
H in terms of matter quantities and we express pt by using the Vlasov 
equation. Integrating by parts and using the constraint equations for pt and 
at lead to (|4^). 


□ 


Let us now define the quantity fj by 


1 


= r/ + — In a 


(45) 


The difference fj{t, 9i) — rift, 62 ) can be estimated by integrating the expres¬ 
sion for po resulting from (|^ and using the energy bound. Since (^) does 


not contain the twist quantities this is exactly as in |A1|. The next step is 


to use the constraint equations (^) and (^) to bound the integral Jgi pdO. 
The net contribution to pt from the twist quantities has the opposite sign 
from that of the other terms. So when obtaining an upper bound for the 
integral it can be discarded and the argument proceeds as in |^]. On the 
other hand the twist terms must be taken into account when obtaining a 
lower bound for the integral. The expression which has to be estimated is 
equal to the twist contribution to the energy up to a factor a. This factor is 
bounded since a is monotone decreasing. Knowing that the difference of p 
at two points 61 and 62 at any given time and its integral at any given time 
are bounded it follows that p{t, 9) < C(t) for some bounded function C(t) 
Remark. In the analysis below C{t) will always denote a uniformly 
bounded function on [fo,t+)- Sometimes we introduce other functions with 
the same property only for the purpose of trying to make some estimates 
become more transparent. 
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Next we show that the boundedness of E{t), together with the constraint 
equation (^4|), lead to a bound on |?7|. For any 9i,92 G *5'^, and t G [tO)i+) 
we get by Holder’s inequality 


\U{t,92) 


Uit,9i)\ = 


rd2 

/ Ue{t,9)d9 

Jdi 


< 



(46) 


The second factor on the right hand side is clearly bounded by 

For the first factor we use the constraint equation (34). This equation can 

be written as 




so that for t G [toT+) 


ft o~‘2‘V 

a~^^‘^{t,9) = / — Pi) + + s^H‘f)ds 

Jto 2 

+a“^/^(to, 6 '). (48) 


Since p > Pi, the integrand is positive and bounded by a multiple of the 
integrand of E{t). Letting C denote the supremum of Q;“^/^(tor) over 
we get 


f a ^/“^dd < f 2sE{s)ds + 27rC < 2E{to){t‘^ — tQ)/2 + 2TrC. 
Je, Jtn 


'01 JtQ 

Hence, for any 6 * 1,02 G 5^ we have 


\U{t,92)-U{t,9i)\ <C{t). 


(49) 


Next Jgi U{t,9)d9 can be estimated using the energy just as in |^| since 
the twist quantities do not play a role in that argument. Knowing that the 
difference of U at any two spatial points and the modulus of its integral 
over the circle are bounded we can conclude that U itself is bounded. These 
arguments also apply to A, since the factor is controlled by the uniform 
bound on U. Bounds on e“^F and e~^P[t also follow from these arguments. 
Indeed, let P = We have from the auxiliary equation 


an 
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expression for Pq in terms of J 2 and we get 


\P{t,02)-P{t,d,)\ = 
r82 


< 2 e 


—i)+2U I 


[ ^ Pg{t,d)dO < [ \e^\J2\de 

J 0\ J 9\ 

[ ' Vae‘^^^-^'^\J2\de < C{t)E{t). 

Je, 


( 50 ) 


Note that fj and U are known to be bounded. Similarily, using the auxiliary 
equation (^) we obtain a bound on 

[ P{t,9)d9 
/si 

in the same spirit as for U. This leads to a bound on P itself which implies 
that e“^r is bounded. A bound on e~8Ht follows if we instead let P = 
Using the auxiliary equations (38-41) we get expressions on 
Pq and Pf. In the expression for Pf we get a term containing At which we 
treat as in the case of bounding U, i,e. we use a Holder argument to bound 
that term in terms of E{t). We also use that e^’^T, A and fj are bounded but 
all ideas have already been used above so we leave them out and conclude 
that e~'8\Ht\ is bounded on [toU+) as well. 

Step 2. (Bounds on Ut^Ue, At, Ag,rjt,fie,at and Q{t).) 

To bound the derivatives of U we use light-cone estimates in a similar way 
as for the contracting direction. However, the matter terms must be treated 
differently and we need to carry out a careful analysis of the characteristic 
system associated with the Vlasov equation. Let us define 


V 

Y 


1 


e , .2 


— o (^i + “^0 ) + + “^0)) 


8P 


,4(7 


y/^UtUg + -^^AtAg, 


and 


X = 


C = 


1 

71 

1 

7^ 


{dt + Vadg) 

{dt - \/adg) 


(51) 

(52) 

(53) 

(54) 


A motivation for the introduction of these quantities is based on similar 


arguments as those given in Step 1, Section 4. For details we refer to [ BCIM ]. 

Remark. We use the same notations, X and Y, as in the contracting 
direction, and below we continue to carry over the notations. The analysis 
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in the respective direction is independent so there should be no risk of con¬ 
fusion. 


By using the evolution equation (pq) , a short computation shows that 
at 


ax+Y) = 


-{X + Y) 


2V2a 

1 / e"^^ \ 

~ X VaXtAg)j 

^ {Ut + ^Ue) ^^4U-2r,^2 ^ 

2 V 2 

+ (A, + + 2ae2’>Sa), 

4\/2 


(55) 


X{X-Y) = 


at 


■(X-Y) 


2 \/2a 

~"v/lt ( ~ ^faUtUg -\- -^^{aA^ — -\/aAtAg)j 

^ {Ut - ^^Ue) ^^4U-2g-^2 ^ 

2 V 2 

4v2 


(56) 


Here k = p — Pi -|- P 2 ~ H 3 . Now we wish to integrate these equations along 
the integral curves of the vector fields x a-nd ^ respectively (let us henceforth 
call these integral curves null curves, since they are null with respect to the 
two-dimensional “base spacetime”). Below we show that the quantity 

Vh(t) := sup X(t, •) + Q^(t), (57) 

6 »gS 1 

is uniformly bounded on [toP+) by deriving the inequality 

W{t)<C+ [ W (s) In W{s)ds. (58) 

Jto 

First we note, exactly as in the contracting direction, that the symmetry 
implies that 

y2(^)eC/(i,e(i)) 
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and 


are conserved. Here V^(t) and 0(t) are solutions to the characteristic 

system associated to the Vlasov equation. From Step 2 we have that U and 
A are uniformly bounded on [to,t+). Hence |H^(t)| and |H^(t)| are both 
uniformly bounded on [to,t+), and since the initial distribution function /o 
has compact support we conclude that 

sup{|u^| + |u^| : 3(s,0) e X 5'^ with f{s,6,v) / 0}, (59) 

is uniformly bounded on [to,t+). Therefore, in order to control Q{t) it is 
sufficient to control 


Q^{t) := sup{|u^| : 3(s, 9) G [to, t] x S^such that/(s, 9, v) 7 ^ 0}. (60) 

Below we introduce the uniformly bounded function 7 (t) to denote estimates 
regarding the variables and . As observed in B there is some can¬ 
cellation to take advantage of in the matter term (p — Pi) which appears in 
the equations for X -\-Y and X — Y above. It is proved there that 

{p-Pi){t,9)<C^{t)\nQ\t) (61) 


In a similar fashion we can estimate P 2 , P 3 and S '23 ■ 

Let us now derive (58). As in Step 2 in Section 4 we integrate the 
equations above for X + Y and X — Y along null paths. For t >tQ integrate 
along the two null paths defined by y and (, starting at (to;^) and add the 
results. In this way the following inequality can be derived: 


sup X{t, ■)< C + C{t) / [1-|-sup A(s, •)] ln(5^(s)](is 


(62) 


'^0 


In doing this it is important to use the fact that and the above 

estimates for matter quantities in terms of InQ^. What is new compared to 
[^] is the occurrence of twist quantities and these can be treated using the 
fact that e“^F and e~^Ht are bounded. 

Let us now derive an estimate for in terms of sup^) X. 

Lemma 2 Let Q^{t) and X{t,9) be as above. Then 

\Q^{t)\^ < C + D{t) [ [{Q^{s))'^ + supX{s,-)]ds, (63) 

Jtr) 9 


where C is a constant and D{t) is a uniformly bounded function on [to,t-i-). 
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Proof. The characteristic equation for associated to the Vlasov equation 
and the constraint equations (^) and (|3^ ) imply that 


(64) 


^{V\s))^ = 2 V\s)-^V\s) = Ti + T2 + n + T4, 
as as 


Here Ti, T 2 and T 3 are expressions given in |A1| and they can be estimated 
just as in the corresponding lemma of that reference. The additional term 
T 4 , which collects together the contributions involving the twist quantities, 
is 


T 4 = -2s{V 


1\2 


g-2r, 


+ s‘^Hf) - 2V^is)[e-^{e^^TV‘^ + sHtV^)] (65) 


It is easily seen that IT 4 I < and this suffices to obtain the 

desired estimate. 


□ 


Combining the estimate for {Q^(t))‘^ in the lemma and the estimate (62) 
for sup 0 V(t,-), we find that W{t) satisfies the estimate (^) and is thus 


uniformly bounded. The constraint equation (32) now immediately shows 
that \r]t\ is bounded by 


2tX + 


te 


-2r; 


(e^^r^ + t^H^) + < C{t)[l + sup V(t, •) + {Q{t)y 

e 


since 

P = [ /du<||/o||oo [ dv <C{Q{t)f. 

4 r 3 J\v\<Q{t) 

Recall that e“^’^(e‘^^r^ + t^H^) is known to be bounded. Thus also rj is 
bounded which implies that both |r| and \Ht\ are bounded. The bound on 
r] also provides a bound on Ifj^l using the constraint equation (^), where 
|Ji| is estimated in terms of Q{t). Analogous arguments show that \at\ is 
uniformly bounded. The uniform bound on X provides bounds on \Ut\ and 
\At\, but to conclude that \Ug\ and lA^I are bounded we have to show that 
a stays uniformly bounded away from zero. Equation (^^ is easily solved, 

Q;(t, 0) = a(to, (66) 


where 

F(t,d) := -2te^^^-^\p - Pi) - 
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which is uniformly bounded from below. Hence jC/el and \Aq\ are bounded 
and Step 2 is complete. 

Step 3. (Bounds on df, ag and r]g.) 

The main goal in this step is to show that the first derivatives of the dis¬ 
tribution function are bounded. In view of the bound on Q{t) we then also 
obtain bounds on the first derivatives of the matter terms p, Jk, Sjk and Pk, 
j, A: = 1,2, 3; j ^ k. Such bounds together with bounds on the 6 derivatives 
of the twist quantities almost immediately lead to bounds on ag and r]g. 
Recall that the solution / can be written in the form 


/(T 0, v) = /o(0(o, t, e, v), V (0, t, 0, v)), 


(67) 


where 0(s, t, 6, v), V{s, t, 9, v) is the solution to the characteristic system 


ds 

dV^js) 

ds 


dV"^ 

ds 

dV^ 

ds 


yO’ 


-ive -Ug + - {rjt - Ut)V^ 

za 

2U^,2t.3 


yo -i{vr-{vr)+ ^y. 


'V V'^ 


-UtV^ - V^Ug- 




-{L-Ut)V^ + V^Ug^ 

s 

q 2U yl 


( 68 ) 


(69) 

(70) 


(71) 


with the property 0(t, t, 9, v) = 9, V{t, t, 9, v) = v. Hence, in order to estab¬ 
lish bounds on the first derivatives of / it is sufficient to bound dQ and dV 
since /o is smooth. Here d denotes the first order derivative with respect to 
t,0 or V. Evolution equations for dQ and dV are provided by the character¬ 
istic system above. However, the right hand sides will contain second order 
derivatives of the field components, but so far we have only obtained bounds 
on the first order derivatives (except for r]g,ag). Yet, certain combinations 
of second order derivatives can be controlled. Behind this observation lies 
a geometrical idea which plays a fundamental role in general relativity. An 
important property of curvature is its control over the relative behaviour of 
nearby geodesics. Let 7 (u, A) be a two-parameter family of geodesics, i.e. 
for each fixed A, the curve u ^ ^{u, A) is a geodesic. Define the variation 
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vector field Y := jx(u,0). This vector field satisfies the geodesic deviation 
equation (or Jacobi equation) (see eg. |HE|) 


Bu^ 


Ryyi', 


(72) 


where B/Bu is the covariant derivative, R the Riemann curvature tensor, 
and 7 ^ := 7 u(n, 0). Now, the Einstein tensor is closely related to the cur¬ 
vature tensor and since the Einstein tensor is proportional to the energy 
momentum tensor which we can control from Step 2, it is meaningful, in 
view of (|7^ ) (with Y = dQ), to look for linear combinations of (90 and 
dV which satisfy an equation with bounded coefficients. More precisely, we 
want to substitute the twice differentiated field components which appear 
by taking the derivative of the characteristic system by using the Einstein 
equations. The geodesic deviation equation has previously played an impor¬ 
tant role in studies of the Einstein-Vlasov system (|RM|, [|Rn|| and [[R13|] ). 


Lemma 3 Let 0(s) = &{s,t,9,v) and V^{s) = V^{s,t,d,v), = 1,2,3 be 

a solution to the charaeteristic system Let d denote dt, do or dy, 

and define 




z^ 


a ^/^(90, 
m.l (VtV^ 

+ V 

yl((y2)2 _ ^^g2C/ y0y2-^3 


+ Ue' 


(y0)2 _ (yl)2 
yly2y3 \ 

+^«(y0)2 _ (yl)2 J 

dV^ + V^Ue de, 

p2U 

dV^ - (V^Ug - V^Ag) de. 

s 


./at (V0)2-(Vi)2 


(73) 


(74) 

(75) 

(76) 


Then there is a matrix A = {aim}, l,m = 0, 1, 2,3, such that 


n := 


satisfies 

‘^ = Aa, (77) 

and the matrix elements aim = CLim{s, 0(s), E^(s)) are all uniformly bounded 
on [to,t+). 
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Sketch of proof. Once the ansatz ( 73)-(pq) has been found this is only a 
lengthy calculation. To illustrate the type of calculations involved we show 
the easiest case, i.e. the term. 


dZ^ 

ds 


fL 

ds 


(dV^ + V^UedO) 


d{^V^) + ^Uede 

ds ds 

+V\Ute + Uee^)dQ + V'^UedA. 

ds ds 


(78) 


Now we use (^) and (^0|) to substitute for dO/ds and dV"^/ds. We find that 
the right hand side equals 


d{-UtV^ - V^Ue-^) + {-UtV^ - 




■)UedQ 


+V\Ute + UeeV^VK)dQ + V^Ue (+ V^d(^) 


VO 


\2^V0 


WO- 


Taking the d derivative of the first term we find that all terms of second 
order derivatives and terms containing ag cancel. Next, since 


- ^/aUed 

we are left with 
dZ"^ 


yo 




1^0 




+ y/^UeV^d TTn = -V^Ue^dV^, 


VO 


— = -{UtV^ + ^Ue^-^)Uede - {Ut + ^Ue^)dV\ 


(79) 


(80) 


Finally we express this in terms of T, Z^, and Z^. Here this is easy and 
we immediately get 

U72 T/l 

— = -{Ut + VaUe^)Z‘^. 

Clearly, the map {dQ,dV^) 1 —> ('k,Z^) is invertible so that this step is 
easy also in the other cases. It follows that the matrix elements a 2 m-, = 

0,1,2,3, are uniformly bounded on [to)^+) (only 022 is nonzero here). The 
computations for the other terms are similar. For the Z^ term we point out 
that the evolution equations (36) and (^) should be invoked and that the 
matrix element uiq contains pg and ag/2a, but they combine and form 
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and that derivatives of e^’^F and e~^Ht appear. The latter terms are easily 
seen to be bounded in view of the auxiliary equations. For example, from 
(1^ we get a bound on dg{e~‘^T), 

A(e-^r) = -2e^+^-^^J2 - e^-’?-^^F^(e-^+^^). 

The right hand side is bounded since U, Ug,r], fj, dgf], and F are all bounded 
as was shown in Step 1 and 2. 

□ 

From the lemma it now immediately follows that |n| is uniformly bounded 
on [to,t+)- Moreover, since the system ([f^)-(|7^) is invertible with uniformly 
bounded coefficients we also have uniform bounds on |cl0| and k = 

1, 2, 3. In view of the discussion at the beginning of this section we see that 
the distribution function / and the matter quantities /?, J^, Sjk and P^, are 
all uniformly bounded. From the constraint equation (^) we now obtain 
a uniform bound on ag by a simple Gronwall argument using as usual the 
identity Finally this yields a uniform bound on rjg since 

Oig 

Vo = Vo-^ 

and a stays uniformly bounded away from zero. 

Step 4- (Bounds on second and higher order derivatives.) 

It is now easy to obtain bounds on second order derivatives on U and A by 
using light cone arguments. We define X and Y by 




Y 


— + ^Ufg) 


,4f/ 




+ C!iA%) 


, 4(7 


— \/aUttUtg + -^^AttAfg, 


(81) 

(82) 


and use the differentiated (with respect to t) evolution equations for U and 
A to obtain equations similar to (|55[) and (^^ . In this case a straightforward 
light cone argument applies since we have control of the differentiated matter 
terms. Ugg and Agg are then uniformly bounded in view of the evolution 


equations 


and (CT). Bounds on second order derivatives on / then 


follows from {7^ by studying the equation for dit. The only thing to notice 
is that fjgg is controlled by (33). It is clear that this reasoning can be 
continued to give uniform bounds on [to,t+) for higher order derivatives as 
well. 


□ 
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5 Proofs of the main theorems 


In this section the analytical and geometrical information obtained in pre¬ 
vious sections is combined to obtain the main results of the paper. 

Proof of Theorem |^. For a spacetime satisfying the hypotheses of the the¬ 
orem we know from Section ^ that a conformal coordinate system can be 
introduced on a neighbourhood of the initial hypersurface Sq corresponding 
to the original data. The direct analogues of the results of Section 5 of 
[ ARR| hold and can be proved by the same arguments. In fact the situa¬ 
tion is slightly simpler since in the case of symmetry there is an obvious 
choice of two Killing vectors while in [ ARR | it was necessary to worry about 
choosing two from a total of three Killing vectors in an appropriate way. Cf. 
also [pCIM ] where this type of argument was introduced for the first time. 
By these results it follows that the region where the solution exists can be 
extended to the past so as to include a Cauchy surface Sa of constant areal 
time. Moreover, either the conformal time coordinate extends to all negative 
values, or R tends to zero as the past boundary of the region covered by 
conformal coordinates is approached. In the first of these cases the region 
covered by the conformal time coordinate includes the entire past of the ini¬ 
tial hypersurface in the maximal Cauchy development, as follows from the 
arguments of Section 5 of |ARR]. Also the past of Sa in that region admits 
a foliation by hypersurfaces of constant R. In that region we can transform 
to areal coordinates. For we can choose a new spatial coordinate 6 so that 
its coordinate lines in Q are orthogonal to that foliation. In the second case 
(where R tends to zero on the boundary of the region covered by conformal 
coordinates) the past of Sa is also covered by areal coordinates. It exhausts 
the past of Sa in the maximal Cauchy development, as will now be shown. 
If the spacetime could be extended to the past there would be a sequence of 
points tending to the boundary of the original spacetime in the extension. 
Along this sequence R would have to tend to zero. However the function R 
which is globally defined on the maximal Cauchy development must tend to 
a non-zero limit along the sequence approaching a point of the extension. 
Thus the existence of an extension leads to a contradiction. It follows that 
in both cases the entire past of Sq in the maximal Cauchy development is 
covered. As a consequence of the results of Section Q the spacetime and the 
areal time coordinate can be extended so that the time coordinate covers 
the interval (Rq, oo). It can then be concluded by the argument at the end 
of Section 5 of lARRj that the entire future of Sq in the maximal Cauchy 
development is covered. 
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□ 


Proof of Theorem [^. The mean curvature of the hypersurfaces of constant 
areal time is 

ilk = — Ut + t~^) (83) 

From the field eqnations it follows that 

7]t-Ut + > tu? -Ut ++ ti^Ut - f (84) 


Hence there are Cauchy surfaces with everywhere negative mean curvature. 
Under these circumstances it was shown by Henkel Q that the initial singu¬ 
larity is a crushing singularity and thus a neighbourhood of it can be foliated 
by CMC hypersurfaces. Given one CMC hypersurface the statement in The¬ 
orem 1 about the range of the CMC time coordinate follows from |R11|. It 
remains to see that the CMC foliation covers the entire future of the initial 
hypersurface. This can be proved by an argument nsed in the case of hyper¬ 
bolic symmetry in [A.RR| which will now be recalled. It is enough to show 
that if p is any point of the spacetime there is a compact CMC hypersurface 
which contains p in its past. Let Si be the Cauchy surface of constant areal 
time passing through p. Equation (^) shows that the mean curvature of 
Si is strictly negative. Hence it has a maximum value Hi < 0. Let S 2 
be the compact CMC hypersurface with mean curvature Hi/2. Then the 
infimum of the mean curvatnre of S 2 is greater than the supremum of the 
mean curvatnre of Si and a standard argument shows that S 2 is strictly 
to the future of Si. Hence p is in the past of S 2 , as required. 


□ 
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